Copyright  2005,  Society  of  Photo-Optical  Instrumentation  Engineers. 

This  paper  was  published  in  SPIE  Proceeding,  Image  and  Signal  Processing  for 
Remote  Sensing  XI,  Volume  5982,  2005  and  is  made  available  as  an  electronic 
reprint  with  permission  of  SPIE.  Single  print  or  electronic  copies  for  personal  use 
only  are  allowed.  Systematic  or  multiple  reproduction,  or  distribution  to  multiple 
locations  through  an  electronic  list  server  or  other  electronic  means,  or 
duplication  of  any  material  in  this  paper  for  a  fee  or  for  commercial  purposes  is 
prohibited.  By  choosing  to  view  or  print  this  document,  you  agree  to  all  the 
provisions  to  the  copyright  law  protecting  it. 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2005 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2005  to  00-00-2005 

4.  TITLE  AND  SUBTITLE 

A  Correlated-k  Based  Ultra-Fast  Radiative  Transfer  (kURT)  Method 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Spectral  Sciences  Inc, 4  Fourth  Avenue,Burlington,MA, 01803 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

10 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


A  Correlated-^  Based  Ultra-Fast  Radiative  Transfer  (ATJRT)  Method 


Prabhat  K.  Acharya*a,  Raphael  Panfilia,  Alexander  Berka,  Steve  M.  Adler- Goldena, 

Marsha  Foxa,  Alan  Wetmoreb  and  Irwin  E.  Alberc 
aSpectral  Sciences,  Inc.,  4  Fourth  Avenue,  Burlington,  MA  01803 
bU.S.  Army  Research  Laboratory,  Adelphi,  MD  20783 
cBoeing  Space  and  Intelligence  Systems,  Seal  Beach,  CA  90740 

ABSTRACT 

A  new  sensor-specific  correlated-A  (c-k)  ultra-fast  radiative  transfer  (RT)  formalism,  /cURT,  has  been  designed  for  fast 
broad-bandpass  scene  simulations  from  UV-visible  to  LWIR  wavelengths.  A  higher  resolution  RT  code  (1  cm"1 
MODTRAN)  has  been  adapted  to  output  1  cm"1  correlated-/:  parameters  for  ozone,  water,  and  the  combined  uniformly 
mixed  species  on  a  pressure-temperature  grid,  which  are  merged  to  form  a  compact  c-k  set  incorporating  the  sensor 
bandpass  response  function.  The  compact  set  is  used  to  compute  bandpass  transmittance  and  radiance  in  near-real  time. 
Scattering  parameters  (molecular  Rayleigh,  clouds  and  aerosols),  blackbody  and  solar  functions  are  cast  as  compact  A- 
dependent  source  terms  and  used  in  the  radiance  computations.  Preliminary  transmittance  results  for  3-5  and  8-12  micron 
bandpasses  and  visible-MWIR  sensors  yield  results  within  2%  of  a  1  cm"1  MODTRAN  calculation  with  a  two-orders-of- 
magnitude  computational  savings.  Applications  include  near-earth  broadband  propagation  and  extinction  calculations  for 
target  detection  and  recognition,  mid-range  tracking,  and  search  and  rescue  operations  from  ground  and  low  altitude 
aircraft. 

Keywords:  AURT,  correlated-A,  instrument  response  function,  sensor-specific  bandpass  /.--distributions,  radiance, 
transmittance,  RT  code,  MODTRAN 


1.  INTRODUCTION 

This  paper  presents  preliminary  results  from  AURT  (Correlated-A-distribution  based  Ultra-fast  Radiative  Transfer),  a 
new  fast  and  accurate  radiative  RT  algorithm  designed  for  scene  generation  and  target  acquisition  by  low  flying  aircrafts 
and  land  vehicles.  Here,  “fast”  is  defined  as  two  or  more  orders  of  magnitude  more  rapid  than  1  cm"1  MODTRAN1 
calculations,  and  “accurate”  is  defined  as  an  absolute  transmission  error  no  greater  than  2%  when  compared  to  bandpass- 
averaged  1  cm  '-resolution  MODTRAN  calculations. 

AURT  utilizes  a  sensor-specific  bandpass  A-distribution  (A  is  the  absorption  coefficient)  on  a  temperature  and  pressure 
grid  to  generate  RT  quantities  (transmittance  and  radiance).  Computation  of  the  A-values  is  time-consuming,  but  is 
performed  off-line  prior  to  use.  As  an  example  of  the  time  savings  gained  by  this  process,  the  1-cm’1  MODTRAN 
correlated-A  (c-k)  procedure  applied  to  the  3-5  |im  window  (1334  cm"1  wide)  requires  22,678  individual  A-distribution 
calculations.  In  contrast,  the  new  algorithm  requires  only  17  A-distribution  calculations  for  the  entire  bandpass,  in 
principle,  a  savings  of  three  orders  of  magnitude  in  computational  load. 

AURT  is  designed  to  perform  rapid  transmittance  and  radiance  calculations  for  a  very  large  number  of  lines-of-sight 
(LOS’s)  and  a  smaller  number  of  model  atmospheres  (e.g.,  as  defined  from  Rawinsonde/radiosonde  data)  with  a  given 
sensor  bandpass.  A  high-resolution  kdb  (A-database)  is  pre-calculated  on  a  grid  of  pressure  and  temperature  values  using 
MODTRAN.  For  each  new  sensor  response  function,  these  data  are  convolved  to  form  a  compact  kdb  ( ckdb )  for  each 
atmospheric  “species”.  For  each  new  model  atmosphere,  the  species  ckdb  A-values  are  (P,  T)-interpolated  and  scaled. 
The  results  are  used  for  the  LOS  RT  computations  by  adjusting  the  layer  A’s  for  path  length. 

The  kdb  files  are  independent  of  any  particular  atmosphere  and  will  be  supplied  with  AURT  software.  They  will  cover  a 
large  enough  spectral  range  to  encompass  all  sensor  response  functions  of  interest.  Once  these  files  are  generated,  AURT 
will  not  require  MODTRAN. 

*For  correspondence  please  contact:  prab@spectral.com;  phone  781  273-4770;  fax  781  270-1161;  www.spectral.com 
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We  have  tested  AURT  with  two  relatively  wide  bandpasses,  one  in  the  3-5  pm  region  and  another  in  the  8-12  pm  region, 
as  well  as  with  a  series  of  substantially  narrower  bands  in  the  visible-MWIR.  A  very  wide  bandpass  means  that  all 
physical  quantities,  not  just  the  absorption  coefficients,  may  vary  substantially  in-band.  The  radiance  algorithm  requires 
modeling  this  spectral  variation  in  quantities  such  as  the  blackbody  function,  solar  irradiance  function,  and  cloud,  aerosol 
and  molecular  (Rayleigh)  scattering  properties.  Among  the  molecular  species,  H20  is  a  special  case,  as  its  self¬ 
broadening  and  self-continuum  components  scale  quadratically  with  concentration  while  its  remaining  line  and 
continuum  components  are  linear  with  concentration.  Finally,  inter-species  correlation,  e.g the  correlation  between  C02 
and  H20  transmittance  in  the  3-5  pm  region,  dictates  that  the  bandpass  may  need  to  be  subdivided  into  two  or  more  sub- 
regions  (not  necessarily  contiguous)  to  achieve  the  desired  accuracy. 

2.  ALGORITHM  DESCRIPTION 

2.1  Overview 

AURT  starts  by  using  a  modified  version  of  MODTRAN  to  output  1  cm"1  A-distributions  for  C02,  03,  H20,  and  the 
remaining  atmospheric  species,  which  are  lumped  together  as  a  single  uniformly  mixed  (UMIX)  species.  A  fixed  C02 
concentration  can  be  included  within  UMIX.  These  A-distributions  are  computed  on  a  grid  of  temperature  and  pressure. 
For  scenarios  dominated  by  the  lower  atmosphere  (<  18  km),  the  temperatures  were  chosen  to  be  305,  280,  255,  230,  205 
and  180  K  (6  values),  and  the  pressures  were  selected  to  be  1.0,  0.82,  0.64,  0.46,  0.28  and  0.1  atm  (6  values).  A  finer  (P, 
T)-grid  will  be  needed  to  accurately  model  03.  The  result  is  36  high-resolution  A-distributions  per  species  for  each 
1  cm"1  bin  of  the  bandpass,  which  can  be  used  for  any  sensor  response  function  within  the  covered  wavenumber  range. 
For  example,  a  database  appropriate  for  the  3-5  and  8-12  pm  regions  would  cover  the  3-12  pm  region  (2500  cm"1  wide), 
and  each  species  kdb  file  would  contain  2500x36x17  =  1,530,000  A-values  assuming  17  A-values  per  1  cm’1  bin. 

Each  high-resolution  kdb  file  is  convolved  with  sensor  response  functions  to  create  compact  kdb  (ckdb)  files  as 
described  below.  Since  this  process  creates  just  17  values  for  the  entire  bandpass,  each  ckdb  file  contains  only  36x17  A- 
values.  Both  the  kdb  and  ckdb  files  are  created  off-line;  only  the  latter  is  used  at  run  time  for  the  RT  calculations. 

Each  species  ckdb  file  may  be  used  to  compute  LOS  transmittance  for  that  species.  Under  the  assumption  of  spectral 
uncorrelation  amongst  the  species,  the  total  in-band  transmittance  is  the  product  of  the  species  transmittances  We  have 
found  that  in  order  for  this  assumption  to  hold,  a  broad  bandpass  band  may  need  to  be  subdivided  into  a  small  number 
(typically  2-3)  of  sub-bands. 

The  geometry  module  computes  layer  species  column  densities  and  average  (P,  T)  for  each  LOS.  The  instrument- 
response  A-values  are  obtained  by  (P,  T)-interpolation  of  data  in  the  ckdb  files  and  then  by  appropriate  scaling  by  species 
concentration  and  layer  path  length.  The  /h  (/=  1 ,  17)  A-values  of  each  layer  of  the  LOS  are  summed  up  and  finally  the 
transmittance  is  computed  by  the  formula; 

j  max~17 

T  =  Y,A8jeXV(-HklayerJUlayer)  (D 

7=1  layer 

where  A gj  is  the  weighting  factor  associated  with  the  /-index,  u  is  the  absorber  amount,  and  the  product  Am  denotes  the 

unitless  quantity  of  optical  depth.  Obtaining  kf  for  the  whole  LOS  by  summing  the  layer  kf  contributions  is  based  on  the 
fundamental  assumption  of  correlation  across  layers.  This  assumption,  which  is  accurate  for  typical  atmospheric 
temperatures  and  pressures,  is  basic  to  c-A  applications. 

Lor  fast  radiance  calculations,  it  is  necessary  to  develop  a  combined  A-distribution  for  all  atmospheric  species,  as  well  as 
A-dependent  parameters  that  enter  into  the  calculation  of  source  functions.  The  procedures  for  doing  this  are  described 
in  Sections  2.4  and  2.5. 

2.  2  MODTRAN  c -k  Implementation 

The  c-A  formalism2'3  was  incorporated4  in  MODTRAN  to  recast  band  model  transmittance  into  a  monochromatic  Beer’s 
Law  form  for  more  accurate  interfacing  with  scattering  codes  such  as  DISORT5.  According  to  the  c-A  formalism,  the 

integral  over  frequency  is  transformed  from  frequency  space  (v)  into  cumulative  probability  space  (g).  In  this 
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transformed  space,  the  function  g  (k )  is  defined  as  the  probability  that  an  absorption  coefficient  ( k )  in  the  defined 
frequency  range  will  be  less  than  or  equal  to  k  .  The  expression  for  transmittance  then  becomes 

T=  — 

Av 

where  Y,j(Agj)=l  and  k(gj)  is  the  average  value  of  k(gj  in  the  Agj  interval. 


_  Anax  -N  _ 

]exp(-kvu)-dv  =  ]exp(-k(g)-u)-dg~  £  exp(-k(gj)-u)-Agj 

j= i 


(2) 


The  ^-distributions  in  MODTRAN  are  generated  from  band  model  parameters  at  the  native  resolutions  of  Av  =  0.1,  1,5 
or  15  cm"1;  the  kURT  formalism  uses  the  1-cm'1  resolution.  In  MODTRAN,  a  database  of  k(g)  as  a  function  of  Doppler 
and  Lorentz  line-width  parameters  is  pre -computed  on  a  discrete  grid  of  A gj,  using  a  Monte-Carlo  method  for  randomly 
distributing  n  number  of  lines  within  A  v;  n,  which  ranges  from  1  to  90  per  1  cm’1,  encompasses  the  range  of  the  band 
model  line-density  parameter  values.  The  A gj  were  initially  chosen  so  that  the  k(g)  distribution  with  the  largest  dynamic 
range  in  k  would  have  each  decade  of  k  covered  by  three  discrete  C values.  This  amounted  to  33  ( k ,  Ag)- pairs  per 
spectral  bin  and  species.  This  number  was  later  reduced  to  17  pairs  with  minimal  loss  of  accuracy. 

During  a  MODTRAN  c-k  execution,  a  set  of  layer  kj  (for  each  spectral  bin)  are  computed  based  on  species-averaged 
line-width,  line-strength  and  line-density  parameters  using  the  values  stored  in  precomputed  c-k  and  band  model 
databases.  At  each  layer  the  individual  species  parameters  are  combined  to  generate  ^-distributions  for  a  single  effective 
pseudo  gas  species.  The  layer  ^-values  are  adjusted  via  a  scaling  factor  to  make  the  layer  transmittances  match  band 
model  theory. 

2.3  kURT  Transmittance  Calculation  -  Single  Species 

We  now  derive  an  expression  for  total  bandpass  transmission  for  a  single  species  beginning  with  the  transmittance  for 
the  /th  single  frequency  bin. 

At  the  outset,  the  notation  needs  to  be  clarified:  Index  i  and  associated  variables  refer  to  individual  1  cm’1  bins  of  the 
bandpass.  Index  j  and  associated  variables  refer  to  c-k  terms  within  an  arbitrary  1  cm"1  bin.  Index  /  and  associated 
variables  refer  to  all  pairs  of  (i,  j)  where  i  runs  over  each  bin  and  j  runs  over  c-k  terms  in  each  bin;  i.e.,  1  runs  over  imaxX 
7max  number  of  c-k  terms  encompassing  the  entire  bandpass.  In  what  follows,  we  introduce  a  new  index  m  which  runs 
over  down-sampled  c-k  terms  for  the  entire  bandpass.  For  ease  of  coding  and  compatibility  with  MODTRAN  the  down- 

sampled  A gm  set  is  chosen  to  be  identical  to  the  Ag-  set,  with  both  having  17  members. 


The  transmittance  for  the  ;th  single  frequency  bin  is  given  by 


J  max 

T,  =  I  exp ■  u)  ■  Ag 


j= i 


(3) 


where  j  is  the  index  of  the  (k,  /l#)-pair  within  the  bin.  Including  an  instrument  response  (filter)  function,  normalized  via 

( (;+ 1/2 )  a  v  A 


\f(y)dv 

'  >  y(i— i/2)av 


1 


(4) 


leads  to  the  bandpass  transmittance  given  by 

7max 

Tf  =X/(  7)=X/,  •  X  A gj  exp( ktju) 

i  i  7=1 

We  define  a  cumulative  instrument- weighted  probability  g  for  the  entire  bandpass  such  that 


Ag,  =  f,  ■  A gj 


(5) 

(6) 


where  i  is  the  index  of  the  (k,  zlg)-pair  within  the  bandpass  (i.e.,  £  runs  from  1  to  imllxX jmax,  and  it  represents  the  total 
number  of  c-k  terms).  The  bandpass  probabilities  sum  to  unity: 
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(8) 


The  equation  for  in-band  transmittance  now  becomes: 

W'imax  _ 

Tf=  Tj  4?/  •  Qxp(ktu) 

t-\ 

where  the  sum  over  t  encompasses  all  combinations  of  (i,  j )  pairs.  The  next  step  is  to  down-sample  the  A:,- values  to  a 
manageable  few  for  the  entire  bandpass.  As  mentioned  above,  for  ease  of  coding  and  compatibility  with  MODTRAN,  the 
down-sampled  set  A gm  set  is  chosen  to  be  identical  to  the  17-member  Agj  set. 


The  process  of  down-sampling  starts  with  creating  a  new  cumulative  probability,  gfk )  for  all  ^-values  of  the  entire 
bandpass.  Sorting  the  A,- values  in  increasing  order,  assigning  to  each  A,  the  weight  Ag:  and  defining  a  new  function 

i 

8i(k,)=  using  the  sorted  set  gives  the  new  cumulative  probability  function  g/(k).  The  cumulative  probability 

Z'=l 

function  is  the  inverse  of  the  newly  sorted  k-g  curve.  Whether  the  A/  are  sorted  or  not,  Eq.  (8)  remains  valid.  Assuming 
that  Eq.  (8)  contains  the  Arvalues  in  increasing  order,  down-sampling  leads  to  the  following  expression  for  7}; 

rnmax  - 

Tf  =  t  Agm-exp(kmu)  (9) 

m= 1 

where  km  is  defined  using  the  newly  created  g/( A )  and  using  the  A gm  discretization  of  the  cumulative  probability: 

-  gm-l  ^Sl<8,„  /gm-1^8l<8,n  I  »I 

k,„=  I  (Ag,-k,)/  2  A g[,  g,  =  IA gr,  gm=£Agm'  (10) 

/=1  /  Z=1  l'=l  m'= 1 


2.4  AURT  Radiance  Algorithm  -  Single  Species 

In  the  simplest  approximation,  optical  quantities  t){v)  other  than  spectral  line  extinction  coefficients  are 
assumed  to  have  little  variation  within  the  bandpass,  so  that  their  effective  values  may  be  represented  by  the  bandpass 


average,  A? ,  weighted  by  the  filter  function: 

-  \f{v)-rty)-dv 
\f{v)dv 


(ll) 


Using  this  procedure  for  including  thermal  quantities  (the  blackbody  function  Bv  )  and  the  solar  irradiance  function  ( Iv ) 

in  the  RT  calculations  provides  good  agreement  with  higher-resolution  algorithms  in  the  visible  and  near  infrared  for 
lower  atmospheric  temperatures.  However,  poor  agreement  is  obtained  with  wide  bandpasses  in  the  mid-wave  infrared 
(MWIR),  where  the  blackbody  function  can  vary  by  an  order  of  magnitude  or  more  and  is  a  significant  contributor  to 
overall  radiance.  Here,  it  is  important  to  incorporate  the  spectral  correlation  between  the  blackbody  function  and  the 
molecular  absorption  coefficients.  A  method  for  doing  this  is  described  below.  This  procedure  can  be  generalized  for 
the  solar  function,  scattering  cross-sections,  or  any  other  frequency-dependent  quantity  needed  for  RT  calculations. 


The  bandpass-integrated  thermal  radiance  at  the  front  of  a  layer  is  given  by 

R  =  J-lBv-(l-Tv)-f{v)-dv  =  TfrtBr(l-fr).Ag  (12) 

Av  i  7=i 

where  Bl Ts  the  blackbody  function  averaged  within  the  zth  bin.  We  can  then  rewrite  Eq.  (12),  using  the  notation  of  the 
previous  section,  as  follows: 

*max '  •/ max  /  \  _ 

R  =  S  ’  Ti  =  exP(“M)  (13) 

/= 1 


Proc.  of  SPIE  Vol.  5982  598217-4 


Note  that  Bf  is  independent  of  subscript  j  and  dependent  only  on  subscript  i.  In  the  above  equation,  we  can  assume  that 

k:  are  sorted  and  we  have  a  new  cumulative  probability  function  g{k).  Down-sampling  the  new  probability  function,  we 
can  rewrite 

R  =  HBm-{l-f;}Agm  (14) 

m 

where  Tm=  exp (—k  u)  ,  km  is  as  previously  defined,  and  Bm  is  defined  as  A -weighted  average  blackbody  function: 


-  Sm-\-Sl<Sm  —  / Sm—\-Sl<Sm 

Bm=  Z  (&grB0  /  X  A?/ 
1=1  /  Z=1 


(15) 


The  resulting  layer  radiance  values  R  from  Eq.  (14)  are  combined  with  their  corresponding  path  transmittances  to  form 
the  total  path  radiance,  as  is  done  in  MODTRAN. 


2.5  Combining  (k,  gj-Values  from  Multiple  Species 

As  noted  earlier,  total  transmittance  may  be  rapidly  calculated  from  the  product  of  transmittances  derived  from  the 
species  c-k  databases.  However,  many  applications,  notably  scattering  computations,  require  a  single  ckdb  set,  for  each 
layer,  that  includes  all  species.  One  may  derive  such  a  ckdb  from  a  convolution  of  the  species  ckdb's,  which  obeys  the 
transmittance  product  (spectral  uncorrelation)  assumption.2  The  convolution  gives  the  joint  probability  distribution  of 
the  sum  of  any  two  absorption  coefficients: 

(kA (m)  +  kB(n )) ■  Sgm  ■  Sgn  (16) 

m  n 

This  leads  to  N~  ( k ,  A g  )-pairs.  To  maintain  the  small  number  of  ^-distribution  computations  at  runtime,  the  convolved 
absorption  coefficients  are  re-sorted  in  ascending  order  and  then  down-sampled.  This  step  must  be  performed  at  run 
time,  as  the  absorption  coefficients  are  a  function  of  particle  density  and  path  length,  which  are  not  known  until  an 
atmosphere  and  LOS  is  defined.  Indeed,  this  is  done  for  each  layer  to  compute  transmittances  (for  radiance 
computations).  For  maximum  speed  advantage,  the  convolution  and  resampling  are  performed  sequentially  for  each 
species;  i.e.,  the  procedure  is  done  on  the  first  two  species,  the  down-sampled  result  is  again  convolved  and  down- 
sampled  with  the  third  species,  and  so  on.  We  have  verified  that  the  results  are  essentially  identical  to  those  from  the 
transmittance  product  method. 

A  complication  of  using  a  combined-species  ckdb  is  the  lack  of  a  general  prescription  for  assigning  ^-dependent  optical 
quantities  dy)  to  the  resulting  (k,  Ag)- pairs.  However,  for  scattering  cross-sections  and  associated  quantities  (i.e,  the 
phase  function  and  single-scattering  albedo)  a  straightforward  solution  exists,  since  the  scattering  originates  almost 
entirely  from  a  single  species,  UMIX.  The  ^-dependence  of  these  quantities  is  defined  by  keeping  track  of  how  the 
UMIX  ckdb  values  enter  into  the  combined-species  ckdb  values.  To  save  computation  time,  this  step  is  pre -calculated 
using  the  ground  level  molecular  densities. 


3.  INITIAL  RESULTS 

3.1  3-5  pm  (2000-3334  cm1)  Single  Species  Transmittance 

Single-species  (UMIX  including  C02,  03  and  H20)  transmittance  calculations  for  the  3-5  pm  bandpass  were  computed 
for  vertical  and  45  degree  off-nadir  paths  between  100  km  and  the  ground.  For  the  1976  US  Standard  Atmosphere,  the 
maximum  absolute  error  is  0.0047  in  H20  transmittance  (for  45°  slant  path),  and  maximum  percent  relative  error  is  less 
than  0.7  %  (the  HJRT  and  MODTRAN,  Version  4,  values  were  0.6932  and  0.6979,  respectively). 

3.2  3-5  pm  Total  Transmittance  Results 

Under  the  uncorrelation  assumption,  the  total  bandpass  transmittance  is  given  by: 

T f  =  T”2°  TfMIX  T°3  (17) 

The  uncorrelation  assumption  has  been  found  to  hold  well  for  bandpasses  or  spectral  bins  that  are  narrow  compared  to 
vibrational  band  widths.  Examples  are  the  spectral  bins  used  in  MODTRAN,  which  range  in  width  from  0.1  cm’1  to  15 
cm"1.  However,  for  £URT  the  uncorrelation  assumption  has  been  proposed  to  apply  to  much  wider  bandpasses,  such  as 
one  covering  the  entire  3-5  pm  region,  a  range  of  over  1000  cm"1.  Using  the  above  formula,  we  computed  the  AURT 
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transmittance  for  the  3-5  pm  bandpass  and  compared  the  results  to  MODTRAN.  The  results,  not  too  surprisingly,  were 
not  within  desired  accuracy;  ATJRT  transmittances  were  higher  by  about  0.03. 

The  reason  for  this  discrepancy  lies  in  the  spectral  variation  of  species  transmittances  shown  in  Figure  1 ,  where  one  finds 
strong  anti-correlation  between  water  and  the  remaining  species,  notably  C02.  We  proposed  and  tested  a  remedy  to  this 
problem,  which  is  to  subdivide  the  3-5  pm  interval  into  two  or  more  spectral  intervals  having  different  transmittance 
relationships.  Very  good  results  were  obtained  with  MODTRAN-calculated  transmittances  using  just  two  disjoint  sub¬ 
bands,  one  corresponding  to  strong  water  absorption  and  the  other  corresponding  to  weak  water  absorption,  as  shown  in 
Figure  1.  For  the  ATJRT  implementation  we  chose  a  slightly  different  and  simpler  division  of  the  3-5  pm  band.  The  first 
region  was  2000-2900  cm'1  (3. 4-5.0  pm),  where  there  is  strong  C02  absorption  and  generally  weak  water  absorption,  and 
the  second  was  2900-3334  cm"1  (3. 0-3. 4  pm),  where  there  is  weak  C02  absorption  and  generally  strong  water  absorption. 

The  total  bandpass  transmittances  for  both  the  vertical  and  slant  paths  are  shown  in  Table  1  and  Figure  2.  The  total 
transmittance  is  calculated  by  weighting,  according  to  their  spectral  width,  the  sub-bandpass  transmittances  of  the  two 
spectral  segments  (weight  =  0.6747  for  2000-2900  cm'1  and  0.3253  for  2900-3334  cm"1).  For  brevity,  we  have  not 
reported  here  the  transmittances  for  the  two  sub-bands  individually.  Interestingly,  we  have  observed  that  the  error  in  total 
transmittance  for  the  tropical  vertical  path  is  smaller  than  that  for  the  2000-2900  cm"1  sub-band.  This  is  not  surprising,  as 
the  sub-band  transmittances  are  unequally  weighted  and  their  errors  are  not  all  in  the  same  direction  or  of  the  same  size. 


O 
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strong  water  region 
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Figure  1 :  Spectral  variation  (transmittance  versus  wavelength  in  microns)  of  MODTRAN  water  and  non-water  species  transmittances 
for  a  vertical  path  from  space  to  ground  using  the  Tropical  atmosphere. 

Table  1:  Total  3-5  pm  (2000-3334  cm"1)  bandpass  transmittances  for  the  vertical  and  slant  path  to  ground. 
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Figure  2:  The  3-5  pm  band  transmittances  using  kURT  and  MODTRAN  for  the  six  MODTRAN  model  atmospheres  for  two  viewing 
geometries.  The  six  atmospheres  are  Sub-Arctic  Winter  (SAW),  Mid-Latitude  Winter  (MLW),  1976  U.S.  Standard  (USS),  Sub-Arctic 
Summer  (SAS),  Mid-Latitude  Summer  (MLS)  and  Tropical  (TROP).  The  two  viewing  geometries  are  a  vertical  path  from  space  to 
ground  and  a  45°  off  nadir  path  from  100  km  to  ground. 

3.3  8-12  (ini  (830-1250  cm1)  H20  Transmittance  Results 

The  8-12  (tm  region  presents  a  special  challenge  for  ATJRT  because  of  the  strong  influence  of  the  water  vapor  self¬ 
continuum,  which  scales  quadratically  with  concentration.  Treating  all  of  the  water  vapor  spectral  components  (Figure 
3)  as  a  single  species  requires  referencing  the  self-continuum  to  a  specific  water  vapor  profile  in  the  kdb.  Some  error  in 
the  self-continuum  amount  is  therefore  introduced  when  this  kdb  is  used  with  a  different  profile.  The  problem  may  be 
avoided  by  treating  the  self-continuum  and  the  remaining  components  (line-t- foreign-continuum)  with  two  separate  k- 
distributions,  each  with  their  own  concentration  scaling.  The  two  ^-distributions  may  be  combined  assuming  either 
spectral  uncorrelation  (i.e.,  treating  them  as  two  independent  species)  or  full  correlation.  There  may  be  error  with  either 
choice,  as  the  actual  degree  of  correlation  is  in  between  these  limits. 


Figure  3:  Water  transmittance,  vertical,  space-to-ground,  components. 
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Figure  4  shows  the  vertical  ground-to-space  transmittance  for  the  two  separate  water  vapor  components  and  for  the  total. 
For  each  component  the  /dJRT  transmittance  compares  extremely  well  against  MODTRAN.  The  product  of  the  ATJRT 
line+foreign  (linear-in-concentration  scaling)  and  self-continuum  (quadratic -in-concentration  scaling)  transmittance 
values  is  in  close  agreement  with  total  FLO  transmittance  calculated  by  MODTRAN;  the  largest  error  is  0.0051  for  the 
Tropical  atmosphere.  Thus,  for  this  bandpass  the  assumption  of  spectral  uncorrelation  between  the  self-continuum  and 
the  other  water  vapor  components  holds  well.  The  alternative,  full  correlation  assumption  yields  errors  that  are  several 
times  larger.  Those  errors  may  however  be  mitigated  by  subdividing  the  spectral  region  into  several  sub-bands. 
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Figure  4:  Modeling  water  transmittance,  vertical,  space-to-ground,  by  treating  the  totality  of  line+tail+foreign  components  separately 
from  the  self-continuum  component.  The  left  panel  shows  the  total  MODTRAN  transmittances  and  AdJRT  transmittances  (product  of 
the  two  /TJRT  components).  The  middle  panel  shows  line+tail+foreign  components,  and  the  third  panel  on  the  right  shows  the  self¬ 
continuum  component.  From  bottom  to  top,  the  six  MODTRAN  atmosphere  models  are:  Tropical,  Mid-Latitude  Summer,  Sub-Arctic 
Summer,  US  Standard,  Mid-Latitude  Winter  and  Sub-Arctic  Winter.  The  order  also  corresponds  to  decreasing  amount  of  water 
column  density;  the  Tropical  atmosphere  has  the  most  water,  and  the  Sub- Arctic  Winter  has  the  least  amount  of  water. 

3.4  Visible-MWIR  Radiance 

Visible -NIR  radiance  calculations  are  in  progress  using  a  modified  version  of  MODTRAN  that  incorporates  the  ckdb 
files  and  ATJRT  elements  to  compute  all  RT  quantities.  In  the  initial  code,  bandpass-averages  are  used  for  the  blackbody, 
solar  and  scattering  terms.  Calculations  to  date  have  been  performed  for  visible  through  MWIR  bandpasss  filters  under 
daytime  conditions  with  model  atmospheres  having  a  wide  range  of  water  vapor  and  visibility  (aerosol  amount). 
Multiple  scattering  calculations  were  performed  using  the  DISORT5  option.  Maximum  error  in  radiance  was  found  to  be 
about  5%,  with  typical  errors  of  2.5%  or  less.  Transmittance  errors  were  consistently  less  than  2%.  As  expected, 
computational  savings  relative  to  1  cm"1  MODTRAN  were  proportional  to  the  number  of  ^-values  employed,  and  hence 
to  the  width  of  the  bandpass.  Detailed  results  will  be  presented  in  a  future  paper. 

4.  CONCLUSIONS 

An  instrument-response -function-specific  correlated-/;  approach,  dubbed  ZTJRT,  has  been  developed  to  rapidly  compute 
bandpass  LOS  transmittances  and  radiances  from  the  visible  to  the  IR.  In  its  full  implementation,  /dJRT  will  be  a 
complete  RT  code  that  includes  scattering  due  to  aerosols  and  clouds.  The  goal  is  a  hundred-  to  thousand-fold  speed-up 
with  2%  accuracy  in  RT  quantities  compared  to  a  1  cm"1  MODRTRAN  run.  Preliminary  results  on  several  wide 
bandpasses  and  narrower  visible-MWIR  instrument  response  functions  indicate  that  both  accuracy  and  speed  goals  are 
within  reach. 
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